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Abstract 

In [1,2] the observed decrease in spectral peak energies of IMFs emitted from hot 
nuclei was interpreted in terms of a breakup density that decreased with increasing 
excitation energy. Subsequently, Raduta et al. [3] performed MMM simulations that 
showed decreasing spectral peaks could be obtained at constant density. In this letter 
we examine this apparent inconsistency. 
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1 Introduction 



In a recent analysis of kinetic energy spectra for intermediate mass fragments 
(Z=2 < IMF < Z~20) emitted in energetic light-ion-induced reactions, it was 
shown that the centroids of the Coulomb peaks systematically shift to lower 
energies and the widths broaden as a function of increasing excitation energy. 
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E*/A [1,2]. This result was interpreted as evidence for a decrease in the aver- 
age breakup density from normal nuclear density po at low excitation energies 
to a constant value of p/po ~ 0.3 for E*A ~ 4 MeV and above. Subsequently, 
Raduta et al. [3] performed MMM (microcanonical multifragmentation model) 
calculations [4] that showed that the observed centroid shifts could qualita- 
tively be explained using a constant density of p/po ~ 0.2. They concluded 
that " ... a decrease in the peak centroids for kinetic energy spectra can be 
observed at low constant density, which is different than published results in 
[1]" . Here we examine this apparent conflict in interpretation. 

First, it is relevant to clarify the method used to fit the experimental spectra 
from which the experimental breakup densities were derived. The parameteri- 
zation of the spectral data was based on the model of Moretto [5] , as adapted 
by Kwiatkowski [6]. The model extends the fission transition state model of 
Nix and Swiatecki [7] to describe all mass divisions, accounting for the evo- 
lution of spectral shapes from Gaussian for fission to Maxwellian for nucleon 
emission. In addition to containing Coulomb barrier, slope(temperature) and 
kinematic parameters, the model has the important advantage of allowing for 
a broadening of the spectral widths with increasing excitation energy to ac- 
count for fluctuations that accompany higher E*/A values, as has long been 
known for fission [8]. 

For present purposes the primary concern is the Coulomb parameter, which 
depends on the average charge separation distance between a given fragment 
and the charge distribution of the remaining system. Using fission fragment 
kinetic-energy systematics [9] as a reference point for the charge separation 
distance at low excitation energies, the parameterization in [1,2] provides good 
fits to IMF spectra at E*/A values below 2 MeV, for which p/po ~ 1. For 
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higher excitation energies, satisfactory fits can only be obtained with a reduced 
Coulomb barrier parameter, which implies a larger separation distance for 
the breakup configuration and hence a lower density. In the excitation-energy 
regime where binary breakup is the dominant de-excitation mechanism, the 
parameterization of [1,2] should provide a satisfactory description of the IMF- 
heavy residue charge separation distance at breakup. For multifragmentation 
events, it is assumed that this model provides a first-order approximation to 
the average Coulomb field experienced by individual IMFs at breakup. 

The top frame of Fig. 1 shows the average Coulomb parameters derived from 
the moving source fits in [1,2]. The fits take into account the significant de- 
crease in source size with increasing E*/A, which coupled with the decreasing 
spectral peak centroids, leads to a near constant behavior at high excitation 
energies. For example, at E*/A = 8 MeV the source is about 75% of the 
target mass due to fast cascade and preequilibrium processes. Shown in the 
bottom frame of Fig. 1 are the extracted experimental densities, which assume 
a spherical breakup geometry. Both the Coulomb parameter and the densi- 
ties are observed to decrease rapidly in the range E*/A = 1 - 4 MeV, after 
which they become constant within experimental error. The upshot is that 
up to E*/A ~ 4 MeV the analysis shows a decreasing density as a function 
of increasing excitation energy - which disagrees with the MMM calculations. 
However, above E*/A ~ 4 MeV the decreasing Coulomb peak centroids lead 
to a nearly constant density, which is the same conclusion the Raduta et al. 
predict. 

In order to place the experimental results and the MMM simulations in per- 
spective, it is useful to compare the IMF multiplicity distributions, as shown in 
Table 1. The experimental distributions for specific multiplicities are shown in 
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Table 1 

Comparison of IMF multiplicities and breakup densities between the experimental 
values of [1,2] and MMM simulations [3]. 



E*/A (MeV) 


1.0 


2.0 


3.4 


4.5 


5.7 


6.8 


7.9 


10.0 


MiMF 


0.06 


0.18 


0.5 


1.7 


2.6 


3.3 


4.0 


4.7 


M/mf(MMM) 


* 


* 


2.8 


4.5 


6.0 


7.6 


8.9 


10.1 


M(exp)/M(MMM) 


* 


* 


5.6 


2.6 


2.3 


2.3 


2.3 


2.1 


< P/Po >(exp) 


1.0 


0.95 


0.6 


0.41 


0.26 


0.32 


0.31 


* 


p//9o(MMM) 


* 


* 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 



Fig. 2 , along with the probabilities for events with Mjmf > 2. The multiplic- 
ity data show that below E*/A ~ 4 MeV, IMF emission is primarily a binary 
process. For higher excitation energies, the multiplicities increase rapidly, sig- 
naling the onset of multifragmentation as the dominant decay mechanism. 
With respect to the present discussion, the salient point is that multifragmen- 
tation does not become a significant decay process until the excitation energy 
has exceeded this threshold and therefore comparison with statistical multi- 
fragmentation models such as MMM should be apphed with caution at lower 
energies. 

The comparison of the absolute values and the multiplicity ratios in Table I 
reveals significant discrepancies between the data and the MMM simulations. 
The IMF multiphcity and the charge distributions are basic observables that 
one expects to be reproduced with any multifragmentation model. Two ob- 
servations are apparent in the table. First, for E*/A below 4 MeV, the MMM 
model increasingly overpredicts the multiplicities as the excitation energy de- 
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creases. Second, above this value the model/data ratio is a constant value of 
about 2.3. By increasing the breakup density in the model, this discrepancy 
could probably be resolved. Thus, it is useful to separate the analysis into two 
excitation energy regimes, above and below E*/A = 4 MeV. 

First, we examine the higher E*/A results. In this regime, the experimental 
results of [1,2] lead to the same conclusion as the MMM model; i.e., the kinetic 
energy centroids shift systematically to lower energies with increasing excita- 
tion energy at constant density. This result has been shown previously with 
the simultaneous SMM [10] model in [11,12], a calculation that also repro- 
duces the experimental charge distributions and the multiplicities. Thus, in 
the regime where multifragmentation is the dominant decay mode, the simu- 
lations presented by Raduta et al, which use a more sophisticated Coulomb 
calculation for multifragment events, reinforce previous conclusions. In order 
for the MMM model to yield more quantitative agreement in this region, its 
parameters need to be modified to reduce the multiplicities by a factor of two 
and to reproduce the charge distributions. In addition, the calculation needs 
to include the significant decrease in source size at higher excitation energy. 

Irrespective of any parameterization, the observed centroid decrease with in- 
creasing E*/A must logically be explained at least in part by a decrease in the 
density, as originally pointed out by Poskanzer [13] . The sequential EES model 
[14] and the metastable mononucleus model [15] have been shown to describe 
the density evolution in this low energy regime. Raduta et al do not address 
this necessary evolution in the breakup density from p/po — 0.2 assumed in 
their model to p/po = 1-0 at low excitation energy. 

In conclusion, for excitation energies where multifragmentation is the domi- 
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nant breakup mechanism, i.e. E*/A > 4 MeV, the MMM simulations of [3] and 
the experimental results [1,2] both show a decrease in the IMF spectral peaks 
with increasing E*/A at constant density. Thus, in this regime, the experi- 
mental results and the calculations are qualitatively self- consistent. For lower 
excitation energies, binary breakup is the dominant decay mode and analysis 
of the spectra with a multifragmentation mechanism cannot account for the 
data. 
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Fig. 1. Top: Dependence of the average moving-source Coulomb parameter < kc > 
as a function of excitation energy. Symbols are as follows: 200 MeV ^He (□); E/A 
= 20-100 MeV ^''N (A); 4.8 GeV ^He (□). Bottom: Average density < p/po > as a 
function of E*/A derived from the kc values in the top panel. 
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Fig. 2. Top: Average number of IMFs as a function of excitation energy for the 8 
GeV/c 7r~ + ^^^Au reaction: observed Nj^/ (circles), corrected for geometry (solid 
triangles) and corrected for both geometry and fragment energy thresholds Mj^y 
(open triangles. Middle: probability for a given number of observed IMFs. Bottom: 
Probability for corrected IMF multiplicity Mi^f > 2 (circles) and Mj^j < 3 (solid 
circles) . 
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